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1 Introduction 



Rare B meson decays, induced by flavor-changing neutral current (FCNC) b — > s transi- 
tions, is one of the most promising research area in high energy physics. Theoretical interest 
to the rare B decays lies in their role as a potential precision testing ground for the stan- 
dard model (SM) at loop level and experimentally these decays will provide quantitative 
information about the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements V t d, V ts and 
Vtb- Besides, these rare decays have the potential of establishing new physics beyond SM, 
such as two Higgs doublet model (2HDM), minimal supersymmetric extension of the SM 
(MSSM), left-right models, etc. 

At present, the main interest is focused on the rare B meson decays, for which the 
SM predicts "large" branching ratios and that can be potentially measurable at working 
B factories and LHC The rare B — > K*£ + £~ (£ = e, /i, r) decays are such ones. At the 
same time this decay constitutes a suitable tool in looking for new physics beyond the SM. 
At quark level the process is described by the b — > s£ + £~ transition. This transition in 
framework of the SM and its various extensions have been extensively investigated [f[f-[]I5| . 

One efficient way in establishing new physics is the measurement of lepton polarization. 
This problem is widely discussed in literature for the b — > s£ + £~ decay [TB|]-|19||. Note that 
all previous studies for the lepton polarization, except the work |Tj|, have been limited 
to SM and its minimal extensions. In [19| the analysis of the r lepton polarization for 
the b — > st + t~ decay was performed in a model independent way. In this work twelve 
(ten local and two nonlocal) four-Fermi operator interactions were introduced in a model 
independent way instead of three independent structures which are present in SM. 

It is well known that the theoretical analysis of the inclusive decays is rather easy but 
their experimental detection is quite difficult. However for exclusive decays the situation is 
contrary to this case, i.e., their experimental study is easy but theoretical investigation is 
difficult. This is due to the fact that the description of the exclusive decays requires form 
factors, i.e., the matrix elements of the effective Hamiltonian between initial B and final 
meson states. This problem is related to to the nonperturbative sector of the QCD and it 
can be solved only by means of a nonperturbative approach. 

These matrix elements have been investigated in framework of different approaches such 
as chiral theory P0| , three point QCD sum rules method PH , relativistic quark model by 
the light-front formalism effective heavy quark theory |23| and light cone QCD sum 
rules 

The aim of the present paper is to perform a comprehensive study the lepton polariza- 
tions and CP-violating asymmetry in the exclusive B — > K*£ + £~ (£ = fi, r) decay in the 
SM and three versions of the 2HDM. Note that this decay in framework of 2HDM (model 
I and model II) were investigated in fl5j (second reference), where Pl and Pt and Afb 
were studied. Here in this work we extend these considerations by studying in model III, 
including normal polarization of r lepton and CP-violating asymmetry, using the current 
limits for the parameters of 2HDM coming from low energy experiments like B°-B° mixing, 
p-parameter analysis and b — > 37 decay. 

The paper is organized as follows. In section 2, starting from a general form of four- 
Fermi interactions we derive model independent expressions for the longitudinal, transversal 
and normal polarizations. In section 3, we apply the above-mentioned general results of 
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the lepton polarizations to SM and to three types of 2HDM (so called models I, II and III). 
A brief summary of our results is presented in this section. 



2 Lepton polarizations 

In this section we present the expressions for the longitudinal, transversal and normal 



polarizations of the r lepton in a model independent way. For this aim we follow |I9 | 
where the matrix element for b — ► st + t~ transition is given in terms of twelve most general 
independent four-Fermi interactions: 
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where Cxx are the coefficient of the four-Fermi interactions. Among them, there are two 
non-local Fermi interactions denoted by Csl and C BR , which correspond to —2Cj in the 
SM. Two (Cll, Clr) of the four vector type interactions (Cll, Clr, C R l and C RR ) are 
also present in the SM (in the forms of Cll = Cg ~ Cio and Clr = Cg + C\o). Clrlr, 
C R llri Clrrl and C R l R l are the scalar type interactions and the last two terms with 
coefficients Ct and Cte are the tensor type interactions. 

For simplicity we will take the mass of the strange quark to be zero and neglect tensor 



type interactions, since it was indicated in |26] that the physical observables are not sensitive 
to the presence of the tensor type interactions. Having established the matrix element 
for the b — > sr + r _ transition, our next problem is calculation of the matrix elements, 
(K* |s 7jU (l ± 75 )6| B), (K* |sz<v<f (1 + 7 5 )6| B) and (K* |s(l ± 75 )6| B), in order to be able 
to calculate the physically measurable quantities at hadronic level. These matrix elements 
can be written in terms of the form factors in the following way: 



(K*{p K *,e)\s^{l± l5 )b\B(p B )) 



m B + m K * 



± ie* (m B + rriK^A^q )^fi(p B + p K *)^{e*q) 



A 2 (q 2 



m B + m K * 



2m 



K* 



q 2 



e*q) A 3 (q 2 ) - A (q 2 ) 



(2) 



(K*{p K *,£)\si<j^q' / {l+ l5 )b\B{p B )) = 

^ upa e w p p K *q°T x {q 2 ) + 2i \e*(m 2 B - m|*) - (p B + p K * )^(e*q) ] T 2 (q 



-2i{e*q) 



% - (Pb+Pk*) 



M 2 

m B — m 



K" 



Uq 2 ) , 



(3) 



2 



where e is the polarization vector of K* meson and q = p B ~Pk* is the momentum transfer. 
In order to ensure finiteness of (2) at q 2 = 0, we assume that A 3 (q 2 = 0) = A (q 2 = 0). To 
calculate the matrix element (K* |s(l ± 7s)6| B), we multiply both sides of Eq. (2) by g M 
and use the equation of motion. Neglecting the mass of the strange quark, we get 



(K*(p K *, e) \s(l ± 75 )6| B(p B )) = —{ T i{e*q){m B + m K ,)A 1 {q 2 ) 

±i{m B -m K »)(e*q)A 2 (q 2 ) ±2im K *(s*q) [A 3 (q 2 ) - A (q 2 )] }. (4) 
Using the equation of motion, the form factor A 3 can be expressed as a linear combination 



of the form factors A% and A 2 (see |2j 
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From Eqs. (1), (2), (3) and (6) we get the following expression for the matrix element of 
the B -> K*£+£- decay 
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where 
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At this point we would like to make the following comment. The main difference from the 
SM case is, we have six different structures (after setting m s = and neglecting tensor 
interactions) in the inclusive channel, while under the same conditions we have only two 
new structures, namely scalar type interactions proportional to and B 5 . On the other 
hand there appears no any new structure in the 2HDM. 

Having established this matrix element, let us now consider the final lepton polarization. 
We define the following three orthogonal unit vectors: 



e L = 



e N 



Pi 

ml 

Vk* x pi 



Wk* x pi | 

e T = e L x e N , (9) 

where f>\ and px* are the three-momenta of the lepton t~ and K* meson, respectively, in 
the center of mass of final leptons. The differential decay rate for any spin direction n of the 
£~ lepton, where n is a unit vector in the t~ rest frame, can be expressed in the following 
form 

dT (n) 1 / dF \ r / „ . ^ _ \ i 

-if = 2 (rf?) J 1 + + p " e " + ■ ] • < 10 » 

where the subscript " 0" corresponds to the unpolarized differential decay rate whose explicit 
form will be presented below. P^, and Pt are recognized as the longitudinal, normal 
and transversal polarizations, respectively. It follows from the definition of unit vectors e*j 
that P T obviously lies in the decay plane whose orientation is determined by the vectors pi 
and p K * and P N is perpendicular to this plane. 

The expression for the unpolarized differential decay rate in Eq. (10) can be written as: 
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The polarizations P L , P/v and P T are denned as: 
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After lengthy calculations we get the following general expressions for the longitudinal, 
transversal and normal polarizations of the t~ lepton (for m s = and neglecting the tensor 
interaction) 
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where A is the expression within the curly parenthesis of the unpolarized differential decay 
rate in Eq. (11). These expressions for the longitudinal, transversal and normal polariza- 
tions are general and model independent (if the tensor interaction is neglected). It follows 
from the expressions of Pt and Pn that they are proportional to the lepton mass and 
therefore they are nonvanishing only for the r lepton. In this work we also analyze the 
CP-violating asymmetry, which is defined as 



A C p(q 2 



where (dT /dq 2 )o is the unpolarized differential decay rate given by Eq. (11) and (dT/dq 2 )o 
is the unpolarized differential decay rate for the antiparticle channel. Note that in SM , CP- 
violating asymmetry is equal to zero (or suppressed very strongly), since all form factors are 
real (see below), Wilson coefficients C7 and C\o are real and only Cg^ contains a strong 
phase. But this strong phase can not lead to CP-violation itself. Using these general 
expressions we can study the sensitivity of the r-lepton polarizations on the new Wilson 
coefficients. Furthermore one can investigate how strongly these polarizations deviate from 
the SM predictions and for which Wilson coefficient this departure is more essential. But 
in the present work we will apply these general results to concrete models, namely to the 
SM and three type of 2HDM, i.e., models I, II (about models I and II, see for example 
27 and III. Note that in models I and II, the flavor changing neutral currents which 

2Sf 



The phenomenological 



appear at tree level are avoided by imposing ad hoc symmetry 

consequence of the 2HDM without this discrete symmetry has been investigated in 
(see also |30|-||39[|). One novel feature of model III is existence of new weak phase which 
appears in Yukawa interaction of fermions with Higgs fields (see below). Existence of this 
new weak phase can lead to sizeable CP violation in B —>■ K*i + i~ decay. Therefore if in 
future experiments sizeable CP violation in the B — * K*£ + £~ decay is discovered, it is an 
unambiguous indication of the existence of new physics beyond SM, since in the SM the 
CP asymmetry suppressed very strongly. 
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Making the following replacements in the expressions given in Eq. (8), the explicit forms 
of Ai, Bi, Ci and Di can be obtained in SM and 2HDM easily. 



1. SM 



Cll = Cg n (m b ) - Cio(m 6 ) , 

C RL = 0, 

C BR = -2C e 7 ff (m b ) , 
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2. 2HDM 
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for example ||0], |TJ) 



10 \ m b) 



Cl" 2HDM (m b )-C™ DM (m b ) 
, 

_ 2C eff2HDM {mb) ) 

C 9 e//2HDM (m fe ) + C 2 ^ M 
o , 

—Cn„ . 



(15) 



2HDM i 



lm W/ 



x ■ 



(7- 5x 



7, 9 and 10) to the leading order are given by (see 
3x 2 ) x 2 (3x-2) 



24(x 



+ 



4(x 



lnx 



+ |Att| 



+ A tf A 



'y(7-5 2 /-8 2 / 2 ) + ^(%-2) 



6b 



72(y - 1)3 12(y 
y(3 - 5y) + 2/(32/ - 2) 



12(y -1)2 6(2/ -i; 



-I) 4 
In 2/ 



(16) 



ri2HDM 
U 9 



+ 



1 



sin 2 9 w 
x 2 (25- 19a;) 
36(s - l) 3 



BH + - — 4 ?r ^ C( m w) 



-3x 4 



sin 2 9w 
f 30x 3 - 



54x 2 + 32x - 8 



+ I A tt I 



1 - 4sm 2 6> 



sin 2 9 



w 



18(x-l) 4 

w xy f 1 1_ 

8 U-l (2/- 



lny 



In a; + 



7 



/47y 2 -7% + 38 3y 3 -6y 3 + 4 n N 
y\ -.no/.. 7^ 7777- 7T77- ln ?/ 



108(y- l) 3 



18(2/ - I) 4 



(17) 



r<2HDM 
°10 



sin 2 9 



w 



(B(m w ) - C(m w )^j 



+ I A t 



sin 2 9 w 8 



+ 



y-1 (y-l) s 



In y 



(18) 



Cqi(^w) = 



+ 
+ 



m^mi 1 



1 x 



m l° |A«| sin 2 9 w 4 v 

+ ^sin 2 a + hcos 2 aj (1 — z) 



j-j |sin 2 a + /icos 2 o;^ f\(x,y) + 



sin 2 2a 
2m 2 H± 



' 2 Ko+m^p) 2 
fc ° 2^ 



/ 3 (y) 



(19) 



where 



C(s) 

f2(x,y) 
My) 



m^mn 1 
|A«| 



rn 



2 ' 



fi(x,y) + 



1 + 



2 2 



m 



w 



f2(x,y) 



X 

y 



h 



m 



hP 



m 



+ 



4(x-l) 4(x-l) 2 



In a; , 



x 



x-Q 3x + 2 , 

+ — 7T77 lnX 



4\^2(x-l) 2(x-l) 2 
x lnx y my 



x - 1 y 
a; lny 



1 ' 
+ 



In ,2 



(z-l)(x-l) 
l-y + y\ny 

(y - 1) 2 ' 



(20) 



(21) 



sin 2 9w = 0.23 is the Weinberg angle, h°, H° and A are two scalar and pseudoscalar Higgs 
fields, respectively. The coefficients X tt and X bb for model I and model II of the 2HDM are: 



X tt = cot (5 , A;,;, = — cot (3 , for model I , 
\u = cot (3 , Ay, = + tan /3 , for model II , 



(22) 



while in model III X tt or X bb is complex, i.e., 

XttXbb = lA^Abfel 
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From Eqs. (16)-(20) we observe that the SM results for the Wilson coefficients Cj (mw), 
Cg M (m w ) amd Cf^(mw) (and correspondingly at \i = m b scale) can all be obtained from 
2HDM results by making the following replacements 



(m w ) 
(m w ) 
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The evolution of the Wilson coefficients from the higher scale \x = mw down to the low 
energy scale \x = m b is described by the renormalization group equation. The coefficients 
C 7 e// (/i), C 9 e// (/i), C 10 (m) at the scale 0(n = m b ) are calculated in |12|, [43| and Cq 1 and 
Cq 2 at the same scale to leading order are calculated in [PJj . The Wilson coefficient Ci is 



2HDM, 
10 

SM, 



m w ). 
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In 



not modified as we move from fi = m w to = m b scale, i.e., Ci (m b ) = C 
order to calculate Cg HDM at m b scale, it is enough to make the replacement C, 
Cg (mw) and then solve the corresponding renormalization group equation. Hence, 
including the NLO QCD corrections, Cg (mi,) can be written as: 
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represents the (9 (a s ) correction from the one gluon exchange in the matrix element of Og, 
while the function g (m c , s) arises from one loop contributions of the four-quark operators 
Oi-Oq, whose form is 



g(yi 
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where yj = Amf/p 2 . The Wilson coefficient Cg receives also long distance contributions, 
which have their origin in the real cc intermediate states, i.e., J/ip, ip', ■ ■ ■■ The J/ip family 
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is introduced by the Breit-Wigner distribution for the resonances through the replacement 

3tt ^ m v T(Vi -> £+£-) 



g (m c , s) -> g (m c , s) - 
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where the phenomenological parameter k = 2.3 is chosen in order to reproduce correctly 
the experimental value of the branching ratio (see for example 



3 Numerical analysis 

In this section we would like to present our numerical results. The main free parameters 
\t, ^bb of the 2HDM are restricted from B — > decay, B°-B° mixing, p parameter and 
neutron electric-dipole moment p6| , that yields \Xbb\ = 50, \X a \ < 0.03. Throughout the 
numerical analysis for the mass of the Higgs bosons we have used m^o = 80 GeV, m#± = 
250 GeV, m A o = 250 GeV and m H o = 150 GeV. 



For the values of the form factors, we have used the results of [25|, where the radiative 
corrections to the leading twist contribution and SU(3) breaking effects are also taken 
into account. The q 2 dependence of the form factors can be represented in terms of three 
parameters as 



F(0) 



1 - a F -C + bj 



2 ' 



where, the values of parameters F(0), ap and 6^ for the B 
Table 1. 



K*£ + £ decay are listed in 





F(0) 


dp 


b F 


Af- 




0.34 ±0.05 


0.60 


-0.023 


At 


^K* 


0.28 ±0.04 


1.18 


0.281 


yB- 


^K* 


0.46 ±0.07 


1.55 


0.575 


T B- 


*K* 


0.19 ±0.03 


1.59 


0.615 


rpB- 

1 2 




0.19 ±0.03 


0.49 


-0.241 


rpB- 

J 3 




0.13 ±0.02 


1.20 


0.098 



Table 1: B meson decay form factors in a three-parameter fit, where the radiative correc- 
tions to the leading twist contribution and SU(3) breaking effects are taken into account 



25] 



In Fig. (1) we present the dependence of the CP-violating asymmetry on q 2 and on the 
weak phase 4> for the B — > K*t + t~ decay in model III, since we have already noted that 
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in SM and in models I and II the CP-violating asymmetry is practically zero. We observe 
that CP asymmetry differs from zero in the region < < 2ti, except at = 0, n and 2n, 
and its value in the region O<0<7r(7r<0< 2ti) is negative (positive). 

Fig .(2) depicts the dependence of the averaged CP asymmetry (A C p) (here and in 
all of the following discussions, by the averaged values of the physical quantities we mean 
integration over q 2 in the region 14 GeV 2 < q 2 < (m# — m^») 2 ) on the weak phase angle 
in model III, taking into account short and long distance contributions. It follows from 
this figure that (A C p) varies in the range (-0.04, 0.04) which is different from zero and 
it definitely is an indication of the existence of new physics beyond SM, since (Acp) is 
practically equal to zero in the SM. 

In Fig. (3), the dependence of Pl on q 2 and the weak phase angle without long distance 
effects in model III is presented. From this figure one can see that, for q 2 > 14 GeV 2 , Pl 
varies in the range (-0.65, -0.8) which is larger than the SM prediction. This is due to 
the fact that the "new" contribution which comes from the charged Higgs boson gives 
constructive interference to the SM results. 

In Fig. (4) we present the averaged longitudinal polarization {Pl) on the weak phase 
angle 0, taking into account short and long distance contributions. For completeness the 
predictions of SM, model I and model II on (P L ) are also presented. It is observed from this 
figure that (Pl) in model III as modulo, is larger than the ones predicted by SM, model I 
and model II. Therefore an observation of |(-Pl)| > 0.65 is another conclusive confirmation 
of the existence of new physics beyond SM. 

In Figs. (5) and (6) we present the dependence of the transversal and normal polar- 
izations of the r lepton on q 2 and on the weak phase angle 0, respectively, without long 
distance effects in model III. The dependence of the averaged transversal and normal polar- 
izations on the weak phase angle, taking into account short and long distance contributions, 
are depicted in Figs. (7) and (8), respectively. For sake of completeness we presented also 
the predictions of SM, model I and model II of the same physical quantity in both figures. 
Fig. (7) clearly depicts that, the prediction of model III on {Pt) as modulo, is approxi- 
mately five times smaller than the ones predicted by SM, model I and model II. However 
the situation is totally different for {Pn), having a range of values in the region (0.10, 0.15) 
in model III, it is approximately two or three times larger than the ones predicted by SM, 
model I and model II. 

Here we would like to make the following remark. It follows from Eq. (13) that Pn 
is defined as the imaginary part of the form factors and of the corresponding Wilson co- 
efficients Cy^ , Cg^ , Cio, Cq 1 and Cq 2 . In SM and C w are real and only has 
imaginary part. On the other side all theoretical methods predict these form factors to 
be real quantities. For this reason, if in future experiments a different value for Pn were 
observed compared to the SM prediction, it is an indication of unambiguous information 
about the existence of the above-mentioned CP-violating phase in theory. 

In conclusion, we have investigated the exclusive B — > K*t + t~ decay in the SM and in 
three different versions of the 2HDM. From the results we have obtained we conclude that 
the combined analysis of the CP-violating asymmetry and r lepton polarization effects are 
very useful tools in looking for new physics beyond SM. 
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Figure captions 



Fig. 1 The dependence of the CP-violating asymmetry A C p of the r lepton on q 2 and on 
the weak phase in model III. 

Fig. 2 The dependence of the averaged CP asymmetry (A C p) of the r lepton on the 
weak phase in model III, taking into account short and long distance contributions. 

Fig. 3 The dependence of the longitudinal polarization Pl of the r lepton on q 2 and on the 
weak phase in model III, taking into account only the short distance contribution in Cg . 

Fig. 4 The dependence of the averaged longitudinal polarization (P L ) of r lepton on 
the weak phase 0, taking into account short and long distance contributions. 

Fig. 5 The dependence of the transversal polarization Pt of the r lepton on q 2 and on the 
weak phase in model III, taking into account only the short distance contribution in Cg^ . 

Fig. 6 The dependence of the normal polarization P N of the r lepton on q 2 and on the 
weak phase in model III, taking into account only the short distance contribution in Cg^. 

Fig. 7 The dependence of the averaged transversal polarization (Pt) of the r lepton 
on the weak phase 0, taking into account short and long distance contributions. 

Fig. 8 The dependence of the averaged normal polarization (Pn) of the r lepton on 
the weak phase 0, taking into account short and long distance contributions. 
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